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Abstract 

We demonstrate that the discrete character of the vibrational spec- 
trum of a small crystal accounts for size dependence of its thermody- 
namic properties and melting temperature. Using a self-consistent 
statistical method [Phys. Rev. B 66, 054302 (2002)] we derive the 
Gibbs free energy of free nanocrystalline plates and calculate the ther- 
modynamic parameters as functions of plate thickness for Cu. 

PACS numbers: 64.70.D-, 65.80. +n 



1 Introduction 

The effect of depression of melting point was first predicted for free nanopar- 
ticles by Pawlow p] from an analysis of their thermodynamic properties. 
Further, it was proved experimentally in transmission electron microscopy 
studies of granular films of Pb, Sn, and Bi [2J. Later, a relationship between 
the melting temperature of small metallic species and their size was estab- 
lished on the basis of results of electron microscopy studies of Sn [3J, Pb and 
In [I], Ag, Cu, Al and Ge [5] thin films. A comprehensive analysis of electron 
diffraction data on melting of Au nanoparticles and their comparison with 
predictions of phenomenological approaches were carried out in Ref. [6]. 
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The electron diffraction methods, providing only the visual picture of 
melting of nanocrystals, allow one to relate the melting temperature to the 
particle size, but they give no information about energetic properties of a 
nanosystem. In order to get more insight into nature of melting of nanocrys- 
tals, these methods are often combined with various nanocalorimetric tech- 
niques, thus providing a possibility of estimation of both the melting tem- 
perature and the latent heat of fusion of nanocrystals [7] [lb]. 

Early theoretical studies of thermodynamics of nanocrystals were focused 
on attempts to explain the strong size dependence of the melting tempera- 
ture. In thermodynamic consideration of an equilibrium between solid and 
liquid phases, an important role was ascribed to the surface of the nanopar- 
ticle [U IU El ESI ES]- Some phenomenological models of melting of free 
nanoparticles employed the effect of surface melting inherent in the most of 
solids. For example, the liquid-layer model [H E7J postulates formation of a 
quasi-liquid layer at the surface of a nanocrystal at temperature below its 
melting point. Among recent theoretical approaches to description of size- 
dependent melting of nanocrystals, we would like to mention the liquid-drop 
model [28J taking into account reduction of the average cohesion energy of 
atoms of a nanoparticle with reduction of its size and using an earlier pro- 
posed relation between the melting temperature and the cohesion energy 
[29]. 

Recently it was reported that thermodynamic properties of nanocrystals, 
such as cohesion energy [UJdS], Debye temperature [TTJ [THl ED] , activation 
energy of diffusion [211 E2] , vacancy formation energy [231 El] etc. display 
also size dependence. These experimental facts point to size influence on 
statistical characteristics of atoms. 

A considerable advance in computational capacity made for the last decade 
has enabled first-principles simulation of dynamical behavior and energetic 
characteristics of systems consisting of more than 10 5 atoms [30] [38], i.e. 
nanocrystals in the so-called mesoscale regime. 

A comprehensive molecular dynamics studies of thermodynamics of nano- 
sized systems was carried out in Ref. [3S] for spherical Cu nanocrystals. The 
obtained results demonstrate that both melting temperature and the latent 
heat of fusion decrease smoothly as particle radius decreases. The analysis 
of energetic and structural properties of nanocrystals revealed existence of 
a rather thin surface layer where average potential energy and root-mean- 
square (rms) displacements of atoms are largely different from that of the 
core region. The rms displacements of atoms in the surface layer increase 
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dramatically at temperature below the bulk melting point of the nanoparti- 
cle T m , indicating melting of this layer. The thickness of the surface layer is 
relatively small and almost independent of temperature. The rms displace- 
ments of bulklike atoms of the nanocrystal rise linearly with temperature, 
and melting of the inner region of the particle occurs when the Lindemann 
criterion is fulfilled. Therefore, the molecular dynamics analysis of thermo- 
dynamic properties of nanocrystals [551 EH] revealed that a nanosized particle 
should be regarded as a structurally and energetically heterogeneous system. 

According to results obtained in Ref. [35], the behavior of the average 
potential energy of atoms of a nanocrystal exhibits essentially nonlinear rise 
near T m . As shown in Refs. [3UI BP] for the bulk solids, such a behavior is 
attributed to evolution of anharmonic instability of the phonon system of the 
crystal which is directly related to the melting transition. 

In this work we study size-dependent modification of the phonon spec- 
trum of crystals and its influence on thermodynamic properties. For the 
sake of simplicity, we restrict ourselves with the case of a nano crystalline 
plate. We show that the vibrational spectrum of such a system is substan- 
tially discrete and size- dependent. The phonon spectrum determines the 
parameters of statistical distribution of atomic displacements of the crystal. 
Thus, the average values of energetic properties of the crystal exhibit size 
dependence. In our opinion, this is an important mechanism responsible for 
depression of the melting point and variation of thermodynamic properties 
of free nanocrystals with decreasing of their size. It is necessary to note that 
influence of discreteness of the spectrum of eigenvalues on energetic charac- 
teristics of many-particle finite systems was considered previously for atomic 
nuclei |JT] and nanoparticles of degenerated semiconductors [12]. In these 
studies it was demonstrated that corrections due to the discreteness of the 
spectrum are proportional to the inverse size. 

2 Free energy of a crystal 

To examine size influence on thermodynamics of nanocrystals, we use a 
self-consistent statistical method [13] developed to compute thermodynamic 
properties of simple solids. The method consists, first, in derivation of a bi- 
nary distribution function of atomic coordinates and, second, in a variational 
procedure of computation of interatomic distance and effective parameters 
of quasi-elastic bonds between atoms of the crystal. It was shown that the 
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phonon spectrum of the crystal determines parameters of the distribution 
function and, therefore, the average values of energetic characteristics of the 
crystal. Therefore, this approach allows one to clarify how size-dependent 
modification of the phonon spectrum of a small crystal affects its thermo- 
dynamic properties. The method has been applied to description of thermal 
characteristics of the rare gas crystals (RGC) [13] [H] and some simple fee 
metals [46J. In the framework of this method we also computed the for- 
mation energy of a vacancy in the RGC as a function of temperature and 
demonstrated that approaching the melting point (which is assumed to be 
directly related to the point of anharmonic instability) is accompanied by 
dramatic reduction of the energy required to create structure defects of the 
lattice [39J EH HZ]- As in our previous studies, we assume the interatomic 
interaction to be pairwise and approximated by the Morse potential 



The parameters A, Rq, and a have been determined so that calculated val- 
ues of interatomic distance, bulk modulus, and sublimation energy at zero 
temperature fitted the corresponding experimental data. These parameters 
are given in Ref. [39] for the RGC and in Ref. [46J for Cu and Ag. 

In this work we restrict ourselves with the high-temperature limit when 
the distribution function of atomic displacements in a simple crystal may be 
represented as a product of one-particle Gaussian functions [39] given by 



where q« is the displacement of an atom from the zth lattice site, r = k^T j A 
is the reduced temperature, c is a dimensionless effective parameter of quasi- 
elastic bond of neighboring atoms, and the coefficient n{r) determines con- 
tribution of the phonon spectrum to the distribution width. At high tem- 
perature it is expanded into a power series as 
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is the de Boer parameter for the Morse potential, M is atomic mass. The 
condition for the applicability of the high-temperature approximation is given 
by inequality r > cA. Coefficients n\ are determined by the phonon spectrum 
of the crystal [13] . For example, 
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where K is reduced wave vector [13], 



M ^ 1/2 



is the reduced frequency of a phonon with wave vector k and branch j, and 
ej X (k) is the projection of a phonon polarization vector on the x axis. The 
integration in Eqs. (jSJ) and (JfJJ) is carried out over the first Brillouin zone. For 
a bulk fee crystal, n — 2, n\ — 5/6. Eqs. (JS])-© were derived in Ref. [13] for 
the bulk crystals, but, according to Ref. [IB], such an expression is applicable 
to any finite system of bound harmonic oscillators. 

In the high-temperature approximation the Gibbs free energy $ of a sim- 
ple fee crystal consisting of iV atoms can be written as [4*3l H6j [39] 



g{T, p, c, b) = — — = - + 3r log — + 
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Here i? is the nearest neighbor distance, b = a(R — Rq) is a reduced lattice 
expansion, z = 12 is the coordination number, p = P/a 3 A is reduced external 
pressure, and v = R 3 / v2 is volume per atom. 

The first term in (jSJ) determines the entropy part of the free energy of 
atomic vibrations in the crystal, the second term represents the average po- 
tential energy of interaction of neighboring atoms. The third term determines 
a contribution to the free energy of the crystal due to the cubic anharmonicity 
of collective atomic vibrations, evaluated in the second order of the pertur- 
bation theory (a 3 w 2.31 for the RGC ESI U3, a 3 » 2.9 for Cu, and a 3 w 2.0 
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for Ag [IB])- The next term takes into account the long-range attraction be- 
tween atoms of a molecular crystal, where e and a are the parameters of the 
Lennard- Jones potential and >c\ = 4.91 for the fee lattice [39J. In the case of 
simple metals, this term should be replaced with an electron gas contribution 
to the kinetic energy of the crystal 1461 . 



3 /3tt 2 \ 2/3 K 1 



where m is the effective electron mass. 



3 Vibrational spectrum of a thin plate 

To determine size contribution to the free energy of a nanocrystal of size h 
{h ^> R), we proceed from the assumption that the size-dependent modifi- 
cation of the vibrational spectrum is most pronounced in its long- wave part, 
with the wave vectors k ~ n/h or kR < 1. Such a case corresponds to 
the elastic vibrations of the medium and is described by wave equations of 
theory of elasticity. The consideration of the size dependence of thermody- 
namic properties of nanocrystals will be carried out for the simplest type 
of nanocrystals, a thin plate. In this case, the dispersion relations have a 
rather simple appearance. A consistent consideration of elastic vibrations 
and propagation of elastic waves in plates was originally performed by Lamb 
[49] in the framework of general theory of elasticity. 

Let us consider a plate of thickness h, with free surfaces parallel to the 
(x, y) plane and the origin taken at the center of the plate. The z axis 
is normal to the free surfaces z = ±h/2. The system is assumed to be of 
macroscopic size in the (x, y) plane, so that vibrations can propagate as plane 
waves in the x and y directions. Note that we do not consider vibrations of 
the plate as a whole. The components of the stress tensor vanish at the free 
surfaces, 



\±h/2 


= 0, 


(T zz\±h/2 


= 0, 


a *>v\±h/2 


= 0. 



(10) 



For such a system, the displacement vector u can be represented as super- 
position of waves of the horizontal polarization, with u x = u z = 0, u y = u, 
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and waves of the vertical polarization with u x 7^ 0, u z 7^ 0, and u y = 0. It 
is important that these two types of waves do not mix at the boundaries, so 
that they can be considered separately (see, e.g. Refs. [4U| [50]). 

Waves of the horizontal polarization are described by ordinary wave equa- 
tions for the displacement u y , 

Sf + Sf+*X = °, (ii) 

where k — tu/ct, u) is an eigenfrequency, q is transverse sound velocity. A 
solution of ( TTTT) in a particular case of a plane wave propagating along the x 
direction is given by 

u y = C cos{(3z - 8) e* x . (12) 

Parameters j3 and 8 are found from both the boundary conditions (I1UI) and 
the symmetry reasons allowing two types of the solutions, a symmetrical one, 
u y (z = h/2) = u y (z = —h/2), and an antisymmetrical one, u y (z = h/2) = 
—u y (z = —h/2). These solutions lead generally to a dispersion relation given 
by 

= c t k n , (13) 

with 

fcn=^ 2 + C 2 +(^) 2 , n = l,2,3,... (14) 

Here £ and ( are, respectively, projections of the wave vector on the x and 
y axes for plane waves propagating in the xy plane. The projection of the 
wave vector on the z axis takes on discrete values, 

In the case of waves of vertical polarization, the displacement is expressed 
in terms of scalar tp and vector ^ potentials, 

u = grady? + rot (15) 

The potential ^ can be chosen so that ^ y = ^ and ^ x = ^ z = 0. Then (p 
and \l/ satisfy the wave equations, 

A(p + k 2 <p = 0, (16) 
A^ + /t 2 * = 0, (17) 
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where k = u/ci, q is longitudinal sound velocity. The solutions of Eqs. (IT5|) - 
(IT7|) are given by 

y> = (7cos(a;s-<5j)e**, (18) 

# = Dcos(/3*-<y t )e**, (19) 

where a and (3 are real values given by 

a = v¥^, (20) 

= v^e. (21) 

In terms of these notations, the boundary conditions (fTOl) are rewritten as 



dip . _ 



= 0, (22) 

z=±h/2 

= 0, (23) 

z=±h/2 



where p = (£ 2 — /3 2 )/(2£). Hence we get the general dispersion relation for 
the acoustic part of the vibrational spectrum of the nanocrystal, 

From the symmetry reasons, the projections of the displacement vectors u x 
and u z should be either symmetrical or antisymmetrical. At the symmetrical 
conditions 

a? = j7i7r, <5° = (2n + 1)tt/2, m, n = 0, 1, 2, . . . 

the dispersion relation (1241) appears as 

a/i p 2 , . 

cot T tan T = --. (25) 

The antisymmetrical condition is specified by 

Sf = (2n + 1)tt/2, 5° = m7r, m, n = 0, 1, 2, . . . 

leading to 

ah 8h p 2 
tan T cot T = --. (26) 




Figure 1: Dispersion curves of longitudinal vibrations of a thin plate for 
n = 1, ... 5. The dash line corresponds to the dispersion relation u = q£ for 
the bulk crystal. 

Note that satisfying the symmetrical condition for u x implies satisfying the 
antisymmetrical condition for u z , and viceversa. 

In Fig. [1] we plotted dispersion curves for the longitudinal characteristic 
vibrations of a thin plate obtained from (|25|) . At a rather small thickness h 
of the plate, the long-wave part (kR < 1) of the vibrational spectrum of the 
nanocrystal differs substantially from that of the bulk crystal, described by 
the linear relation between the frequency and the wave vector, Uj = Cj£ (j = 
l,t). The spectrum of the thin plate splits to separate vibration branches. 
For each branch, there is an integer number of wavelengths going across the 
plate thickness. The dispersion curves are well approximated by parabolas 
and can be analytically represented as 

^(e,0 = < 6 + «n C/ 2 (e 2 + C 2 ). (27) 

Here is the minimal value of the frequency of the nth vibration branch, 
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a n « 1. For each of both symmetrical (I25p and antisymmetrical (I26p cases 
we have two types of solutions for uj n ^ corresponding to longitudinal and 
transverse standing waves with the frequencies generally written as 

,• 7rn . 

<6~X Cj '' j = Z '*' (28) 

As in the case of waves of the horizontal polarization (fl4|) . the wave vector 
for the both longitudinal and transverse waves of the vertical polarizations 
is given by Eq. ( fl4l) . too. 

A similarity of the expression for (u 3 n b ) 2 to the formula for the energy of 
a quantum particle in a one-dimensional potential well [5 1] is explained by 
the wave nature of both quantum particles and vibrational modes, exhibiting 
discreteness of the eigenvalues in the case of a small size of the system. 

For the bulk fee crystals, the first Brillouin zone has the form of a 
truncated octahedron and may be well approximated by a sphere of radius 
k° s = (67r 2 /f) 1 / 3 . The radius is determined by the requirement that the 
number of independent values of the wave vector falling within the sphere is 
equal to the number of atoms of the simple lattice. 

In the case of nanocrystals, evaluation of the radius k s of the spherical 
Brillouin zone requires taking into account the discrete character of k z in 
replacing summation over k by integration. This can be done using the 
Euler-Maclaurin formula 



(29) 



kmax rkmax 1 

f( k ) = / /(*) dt + 2 + /(AWe)] + 



k 

Then we get 

b — b° I i j , 

4:(k° s R) h J ' 

Presence of a size- dependent factor in (1301) results in appearance of similar 
corrections to the coefficients ni in 

n,(/0=n?(l + 7,^J, (31) 

where is the bulk value of ni, the first two coefficients 7; are 70 = 0.68 and 
71 = 0.74. As will be shown below, the size dependence of the coefficients n\ 
leads to a number of physical effects that should be observed in nanocrystals. 
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It should be noted that the above procedure of calculation of the coeffi- 
cients 7; seems to give overestimated values, because Eq. (I3"UI) for the quasi- 
discrete wave vector is valid, strictly speaking, only in the range kR < 1, 
sensible to the particle size. Near the boundaries of the Brillouin zone, the 
values of k are determined by the structure of the crystal lattice and are 
slightly dependent of the crystal size. However, the structure of Eq. f[3"0l 
should remain, on the whole, the same, so that 7/ may be considered as 
parameters of the theory. 



4 Thermodynamic properties of nanocrystals 

The equilibrium values of the quasi-elastic bond parameter c and reduced 
lattice expansion bo are determined from minimization of the Gibbs function 
(JS]) with respect to c and b, 



dg_ 

dc 



T,p,b 



* 

' db 



0. 



(32) 



T,p,C 



The condition for the equilibrium values of the variational parameters to 
exist is that 



D = det 



g'L 



9cb 
9l 



> 0. 



(33) 



The sufficient condition D = defines the point of anharmonic instability 
of the crystal where the minimum of g(r,p,c,b) with respect to c and b 
disappears. 

In the case of nanocrystals, the Gibbs free energy depends, besides tem- 
perature t and pressure p, on the crystal size h. As a result, a number of 
thermodynamic characteristics of the crystal display size dependence. For 
example, the Debye temperature is directly proportional to the equilibrium 
value of the parameter cq(t, p, h) and is written as [HJ H5] 
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D 



(6v / 27r 2 ) 1/3 (/t ik )^Ac (r,p, h)/k 



(34) 



where coefficients = io^j (kR) account for different polarizations of the 
acoustic waves, and (ftjk) ~ 0.67 for the fee crystals. The volume and tem- 
perature dependence of Qd determines the Griineisen parameter, 



7g 



(d\ne D /d\n V) T 
- (01ne D /01nT), 



(35) 
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Figure 2: Size dependence of the interatomic distance in Cu nanocrystals at 
different temperatures. 

The equilibrium value of bo(r,p,h) governs the interatomic distance, the 
coefficient of thermal expansion, and the bulk modulus of the nanocrystal. 
All these properties would depend on the crystal size. 

If the melting temperature T m is assumed to be identical to the temper- 
ature of anharmonic instability of the crystal [39l HO] , then it can be shown 
in the framework of the present statistical model that T m is related to the 
factor n(r) as 

k B T m 1 v 

r m = — ; 36 

Hence it follows that even a slight increasing of n(r) accompanying reduction 
of the crystal size results in appreciable depression of the melting temperature 
of the nanocrystal. Since R/h 1, it follows from (|36|) that 
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where is the melting temperature of the bulk material, and j3 = 97 -R. 
This is the well-known relationship between the relative melting temperature 
of a plate nanocrystal and its inverse size (see, e.g., [22] )■ The calculated 
value P = 1.59 nm for Cu is somewhat higher than that extracted from 
experimental data [28] (/3 ~ 0.9 — 1.2 nm). 

In Fig. [2] we plotted the interatomic distance in Cu nanocrystals versus 
particle size calculated from (j3"2"j) at three values of temperature. Hereafter 
all theoretical results are presented for zero external pressure. A nanocrystal 
of size h ~ 1.3 nm melts at T = 960 K, which explains some nonlinear rise 
of R(h) in the top curve. These results comply with the assumption [5H [55] 
that the size-induced lattice contraction observed in spherical nanoparticles 
is due to capillary pressure only. 




Figure 3: Size dependence of the effective quasi-elastic bond parameter 
co(r,h) of Cu nanocrystals at different temperatures. 

While the interatomic distances in the considered nanocrystals show no 
marked size dependence, the effective parameter c (r,p,h) of the quasi- 
elastic bond varies substantially with h. Such dependence calculated for Cu 
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nanocrystals is illustrated with Fig. [3] which demonstrates that c decreases 
nonlinearly as h decreases even at T T m (h). 

Size dependence of the Debye temperature Qd of small crystals was ob- 
served in a number of experimental studies [T71 dH [20]. In the framework of 
the present approach, 0£> is proportional to the quasi-elastic bond parameter 
c (154"|) . The Qo(h) dependence presented in Fig. 0] agrees qualitatively with 
that obtained experimentally, e. g. for gold nanoparticles [T9"] . 




Figure 4: Size dependence of the Debye temperature of Cu nanocrystals at 
different temperatures. 

In many bulk solids, the melting transition is preceded by anomalous be- 
havior of temperature derivatives of thermodynamic functions (isobaric heat 
capacity, thermal expansion coefficient, Griineisen parameter) due to evolu- 
tion of the anharmonic instability in the phonon subsystem of the crystal 
as T — ► T m [39l HQ] (premelting phenomena). As seen from Fig. [5] for the 
isobaric heat capacity Cp, similar behavior should be also observed when 
the melting point is approached via reduction of the crystal size at constant 
temperature (the top curve). 
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Figure 5: Size dependence of the isobaric heat capacity of Cu nanocrystals 
at different temperatures. 



In Fig. [6] we show the melting temperature of a free Cu thin plate as a 
function of its thickness calculated from Eqs. (1521) and the condition D — 0, 
compared with that obtained from Eq. (|36j) . 

At present, the only generally recognized melting criterion is the empirical 
Lindemann rule suggesting that a solid melts when the rms displacement of 
atoms reaches a characteristic fraction of the interatomic distance. In terms 
of notations accepted in this work, the Lindemann ratio 5 is given by 



3r 



4cg(r)(ai2)V 



(38) 



Here g t ~ 0.77 is the correlation smearing of the distribution width of an 
atom at a lattice site pl3] . 

The size dependence of the Lindemann ratio is presented in Fig. [7J A 
dramatic increase of the rms displacement of atoms in the premelting region is 
caused by nonlinear reduction of the quasi-elastic bond parameter cq(t,p, h) 
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Figure 6: Melting temperature of a free Cu thin plate versus its thickness 
calculated by minimization of the Gibbs free energy dHJ) (solid curve) and 
from Eq. (|36|) (dotted curve). 

in the vicinity of the instability point (Fig. [3D- The maximal value 8 ~ 0.1 
depends only slightly on h. Thus, as for the bulk solids, the premelting 
region is characterized by a nonlinear rise of some solid-state parameters of 
nanocrystals (isobaric heat capacity, the coefficient of thermal expansion, rms 
displacement etc.). It was also shown for the bulk case that the formation 
energy of the structural lattice defects exhibit a sharp drop just before T m . 
All these premelting effects promote transition to a structurally disordered 
phase, i.e. melting. 

Fig. [HHTUl illustrate nonlinear behavior of thermodynamic properties of 
nanocrystals when the instability point is approached via increasing temper- 
ature at constant size. We see that the premelting range shifts as the crystal 
size decreases. 

If two nanocrystals of the same thickness h are brought into tight me- 
chanical contact, the vibrational modes with wavelengths larger than the 



16 



0,105- 
0,100- 
0,095- 
0,090- 
0,085- 

— i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i 
2 4 6 8 10 12 14 16 18 

h (nm) 

Figure 7: The relative rms displacement of atoms (the Lindemann ratio) of 
Cu thin plate versus its thickness at different temperatures. 

characteristic thickness of the contact become collectivized, and their spec- 
trum becomes, roughly speaking, equal to the spectrum of a nanocrystal 
with thickness 2h [53]. Breaking the contact, i. e. separating the two plates, 
results in the reverse rearrangement of their vibrational spectra and is accom- 
panied by a change Ash in the energy per atom e = g — rdg/dr. Obviously, 
Ash = —Ae 2 h- If the nanoparticles are isolated, the energy change would 
lead to a change in the temperature of the nanocrystals, AT = Ash/cy, 
where cy is isochoric heat capacity per atom of the nanocrystal. In Fig. [TTJ 
we plot a contact temperature change AT of two identical Cu thin plates 
versus initial temperature calculated for two values of the plate thickness, 
h = 5.1 nm and h = 1.26 nm. As follows from the calculations, the tem- 
perature change is maximal if the initial system of separated nanoparticles 
is close to the melting point. 
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Figure 8: The temperature dependence of the isobaric heat capacity of Cu 
nanocrystals: (a) bulk crystal; (b) h = 10.2 nm; (c) h = 5.1 nm; (d) h = 
2.55 nm; (e) h = 1.26 nm. 



5 Discussion 

So far there is a large amount of experimental data concerning size influence 
on the bulk properties of nanoparticles, such as cohesion energy [T7], [18] , De- 
bye temperature (TOES [20], and activation energy of diffusion [2UI22]- One 
of the most pronounced manifestations of change of thermodynamic proper- 
ties observed in nanosized crystalline systems is the effect of reduction of the 
melting temperature of free nanocrystals [2]-[16j. So the principal problem 
of statistical description of thermodynamics of nanocrystals is elucidation of 
the mechanism of size effect on statistical characteristics of atoms in such 
systems. Relying on the results of molecular dynamics simulations of ther- 
modynamic properties and melting of nanocrystals [35l [36] , it is natural to 
conclude that a thin (relative to the particle size) surface layer has only a 
minor influence upon the bulk properties of the particles; moreover, in the 
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Figure 9: The temperature dependence of the isochoric heat capacity of 
Cu nanocrystals: (a) bulk crystal; (b) h = 10.2 nm; (c) h = 5.1 nm; (d) 
h = 2.55 nm; (e) h = 1.26 nm. 

case of spherical particles an additional capillary pressure should contribute 
to increasing of both the melting temperature and the Debye temperature. 

It is shown in this work that an important size- dependent factor govern- 
ing thermodynamic behavior of a nanocrystal is discreteness of its vibrational 
spectrum. It leads to increasing of the parameter n(r) of the statistical dis- 
tribution function of atomic coordinates (j2J) as crystal's size decreases. As 
a consequence, the average value of the interaction energy of atoms in the 
nanocrystal changes with its size. It is necessary to emphasize that such 
influence on a crystalline system mediated by direct change of the statistical 
distribution function is inherent only in nanosystems. In this connection, it is 
worth noting a specific character of thermodynamic response of a nanocrys- 
tal to variation of its size. As the crystal's size decreases, the parameter 
c (t, p, h) of quasi-elastic bond of its atoms decreases (Fig. ED, as well as the 
Debye temperature (Fig. 0]), while interatomic distance remains nearly con- 
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Figure 10: The temperature dependence of the Griineisen parameter of Cu 
nanocrystals: (a) bulk crystal; (b) h = 10.2 nm; (c) h — 5.1 nm; (d) h = 
2.55 nm; (e) h = 1.26 nm. 

stant (Figs. [2]). The melting temperature (I36I) of the nanocrystal is markedly 
decreased (Figs. ED- Independently of the crystal size, the melting transition 
occurs when the Lindemann criterion is satisfied (Fig. [7]). Along with a size- 
dependent shift of the melting temperature, there is also a corresponding 
shift of the premelting range where thermodynamic properties (isobaric heat 
capacity, coefficient of thermal expansion etc.) display nonlinear behavior. 
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Figure 11: A temperature change at mechanical contact of two identical Cu 
thin plates of thickness h versus initial temperature. 
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